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Abstract 

This paper studies the cooperative global robust stabilization problem for a class of nonlinear multi-agent systems. The problem 
is motivated from the study of the cooperative global robust output regulation problem for the class of nonlinear multi-agent 
systems in normal form with unity relative degree which was studied recently under the conditions that the switching network 
is undirected and some nonlinear functions satisfy certain growth condition. We first solve the stabilization problem by using 
the multiple Lyapunov functions approach and the average dwell time method. Then, we apply this result to the cooperative 
global robust output regulation problem for the class of nonlinear systems in normal form with unity relative degree under 
directed switching network, and have removed the conditions that the switching network is undirected and some nonlinear 
functions satisfy certain growth condition. 
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1 Introduction 

Consider a class of cascade-connected multi-agent non¬ 
linear systems as follows: 

Z^ = Fi{Z^,ei,d{t)) 

Cj. — Gii^Zi, €-i, -j- biUi, % — 1, • * • , If 

where {Zi, et) G x R is the state, Mi G R is the in¬ 
put, Ci is the output, hi is an unknown real number, and 
d : [to, oo)—C R””^ with D a known compact subset 
represents external disturbance and/or parameter vari¬ 
ations. It is assumed that < \bi\ < bM for some 
known positive real numbers bm and bM-, and the func¬ 
tions Fi : R"* X R X R"'^—>-R"'*, and Gi : R"‘ x R x R”'^—>-R 
are both sufficiently smooth and satisfy Fi (0,0, d{t)) = 0 
and Gi(0,0, d(t)) = 0 for all d(t) G R"''. 

To describe our control law, let cr : [0, oo) -T V = 
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{1, 2, • • • , no} for some integer no > 0, be a piecewise 
constant switching signal, and Save[Td, -^o] be the set of 
all signals possessing the property of average dwell-time 
Td with chatter bound Nq [8,9], FIp = [h^-] G R^^^, 
p = 1, ■■ ■ , no, be some Ai matriccHI. Dehne a piecewise 
switching matrix Jda-(t) C such that, over each in¬ 
terval [ti, ti-i-i), = Hp for some integer 1 < p < no- 

Denote the elements of by hij{t), = 
and define the virtual output of (1) as 

N 

^va{t)i ^ ^ ^ij ^ 1, * ' ' , ^• (2) 

i=i 

Then, we describe our control law as follows: 

Ui — ^2 (c^cr(t)i), t — 1, Af (3) 

where the functions fci(-), i = 1, • • • , N, are sufficiently 
smooth vanishing at the origin. Such a control law is 
called a distributed switched output feedback control 
law, since we can only use e„o-(t)i instead of for feed¬ 
back control due to the communication constraints, 
which will be further elaborated in Section 4. 


^ A matrix M G is called an M matrix if all of its non¬ 

diagonal elements are non-positive and all of its eigenvalues 
have positive real parts. 
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Problem 1 Given the multi-agent system (1), a set of 
Ai matrices Hp G , p = 1, ■ ■ ■ , uq, and some com¬ 

pact subset D C M"'' with 0 € D, find Td, Nq, and a con¬ 
trol law of the form (3) such that, for any d(t) S D, and 
any a € Savei^d, ^o]: ihe equilibrium point of the closed 
loop system composed of (1) and (3) at the origin is glob¬ 
ally asymptotically stable. 

The above cooperative global robust stabilization prob¬ 
lem is of interest on its own and has not been studied 
before, since (3) is a switched control law which results 
in a switched closed loop system. On the other hand, it 
is motivated from the study of the cooperative global ro¬ 
bust output regulation problem for a class of nonlinear 
multi-agent systems with switching network. In fact, we 
will show that our main result is applicable to the coop¬ 
erative global robust output regulation problem for the 
unity relative degree nonlinear multi-agent systems with 
switching network in Section 4. It is also noted that the 
result in 4 includes some existing results as special cases. 
If no = 1, then reduces to a constant matrix. For 
this special case, the above problem has been studied in 
[ 2 ] under the assumption that Hi is a symmetric A4 ma¬ 
trix. More recently, for the case where no is any positive 
integer, the above problem was studied in [9] under the 
assumption that Hp, p = 1 , • • • , no, are all symmetric 
A4 matrices and the nonlinear function Gi satisfy cer¬ 
tain growth condition |3]. In this paper, we will further 
study the above problem for the more general case where 
the nonlinear function Gi does not satisfy any growth 
condition, and Hp, p = 1, ■ ■ ■ , no, are any Ai matrices, 
i.e., do not need to be symmetric. As an application of 
our main result, we will obtain the solution of the same 
problem studied in [9] without assuming that the non¬ 
linear functions Gi satisfy any growth condition, and the 
communication graph of the multi-agent system is undi¬ 
rected for alH > 0 . 

Without assuming the growth condition, and the undi¬ 
rectedness of the communication graph, the problem in 
this paper is technically much more challenging than the 
one in [9]. To overcome these difficulties, we need to de¬ 
velop a changing supply pair technique for exponentially 
input-to-state stable (exp-ISS) nonlinear systems to re¬ 
move the growth condition of the nonlinear function Gi, 
and we need to employ a non-quadratic function for the 
closed-loop system in Section 3 to handle the directed 
communication graph. 

Notation. For any column vectors Oi, z = I,..., s, denote 
col(ai,..., Os) = [af, . A function p{-) is called 

SAf function if it is a smooth non-decreasing function 
p : [ 0 , oo)—>-[ 0 , oo) satisfying p{s) > 0 for all s > 0 . 
d(s) = 0(a(s)) as s ^ 0+means lim sup < cxd. The 

S-J-0+ 

notation Amin (A) denotes the minimum eigenvalue of a 
symmetric real matrix A. 

^ See Assumption 8 for the definition of growth condition. 


2 Preliminaries 

In this section, we will establish a technical lemma. Con¬ 
sider a general nonlinear system: 

x = f{x,u) (4) 

where x € K" is the state, u G M"* is the input, / : M" x 
is locally Lipschitz and f(0, 0) = 0. It is known 
from [II] that system (4) is said to be input-to-state 
stable (ISS) if there exists a G^ function V : i?”—>-[0, oo) 
such that, for all x,u, 

aidkll) < V{x) < Q! 2 (||a;||) 

^f{x,u) <-a{\\x\\)-G-f{\\u\\) 

for some class /Coo functions oi(-)! o:(-) and some 

class/C function 7 (•). The pair offunctions ( 7 , 0 ) is called 
a supply pair for system (4), and V (a:) is called an ISS 
Lyapunov function of (4). Moreover, if system (4) is ISS, 
then, for any class /Coo function q;(-) satisfyiiig q;(s) = 
0{a{s)) as s —> 0+, there exists G^ function V{x) such 
that, for all x,u, 

didkll) < V{x) < d 2 d|a;||) 

dV (6) 

-^f{x,u) < -Q:d|a;||) -b7d|u||) 

for some class /Coo functions di(-)i and some class 
/C function 7 (-) [II]. This result is called the changing- 
supply pair technique which plays a key role in finding 
a suitable Lyapunov function for a nonlinear system to 
conclude the asymptotic stability of its origin. However, 
as will be explained in Section III, this version of the 
changing supply pair technique is not adequate for han¬ 
dling the stability of switched systems. We need to fur¬ 
ther establish a lemma for the following class of nonlin¬ 
ear systems: 

X = f{x,u,d[t)) (7) 

where x G K" is the state, u G K™ is the input, d : 
[to,cxD)—i>D C with D some non-empty set, repre¬ 

sents external unpredictable disturbance and/or internal 
parameter variation, / : K" x R™ x D— is piecewise 
continuous in d{t) and locally Lipschitz, d{t) is piecewise 
continuous in t and /(0,0, d{t)) = 0 for any d{t) G R"”^. 

Lemma 2.1 Suppose that there exists a G^ function 
V{x, t) such that, for any d(t) G D, 

aidlxll) < H(x,t) < 02(11x11), Vx (8) 

V{x,t) <—XV{x,t)-\-Piu), y x,u (9) 
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where q:i(-) and a 2 {-) are some known class ICoo func¬ 
tions, X is some known positive real number and /3(-) is 
some known smooth positive definite function. Then, 

(i) For any function p G SM, the following function 

rV{x,t) 

V{x,t)= / p{s)ds, (10) 

^0 

is an ISS Lyapunov function in the sense that 

aidkll) < < Q; 2 (||a:||), V a; (11) 


(1) ^ In this case, 

p{V(x,t)){ -\V{x,t) +/3(u)) 
<p{V{x,t)){ - \V{x,t) + ^V{x,t)) 

= - ^p{V{x,t))V{x,t) 

< - [x, t))V (x, t) + p{jl3{u))/3{u). 

(2) ^V[x, t) < j3{u)\ In this case, we have 


V{x,t) <—^V{x,f)-\-j3{u), V X, u, V(i(t) € D (12) 

for some class Koo functions ai(-) and a 2 {-), some 
smooth positive definite function (d{-). 


(ii) For any class ICoo function a(-) satisfying a(s) = 
0(Qfi(s)) as s ^ 0+, there exists some function p G SAf 
such that, for any d(t) G D, V(x,f) satisfies (11), and 
the following 


V{x,t) < - {XV{x,t) + ca{\\x\\)) + I3(u), M x,u (13) 


for any positive real number 0 < A < — and positive real 

number c = -A. 

2 

Proof: Part (i): It is easy to see that 


V{x,t) < p{V{x,t))V{x,t), (14) 


p{V{x,t)) < p{jfi{u)) 
which implies that 

p{V (x, t)){-XV (x, t) + P{u)) 

< - Xp(y (x, t))V (x, t) + p{jP{u))/3{u) 

< - ^p{V (x, t))V (x, t) + p{jl3{u))/3{u). 
Choose a smooth positive definite function fi{-) such that 

fi{u) > pijfiiumu) (17) 


From (14) , (16), and (17), we have 

V (x, t) <- ^p(y (x, t))V (x, t) + fi{u) 


(18) 


and equations (8) and (10) imply that 




p{s)ds < V{x,t) 


< 


ra2{\\x\\) 

/ p{.s)ds < a 2 i\\x\\), Vx 

Jo 


(15) 


Part (ii): Since q;(s) = 0(ai(s)) as s —>■ 0+, by Lemma 
2 of [11], it is always possible to find a function p G SAf 
such that, for all x, 

p(ai(||x||))ai(||x||) > d(||x||). (19) 


Thus 


for some class ICoo functions ai(-) and a 2 (-)- 


p{V{x,t))V{x,t)>p{ai{\\x\\))ai{\\x\\) > a(||x||). (20) 


We now show that along the trajectory of x — Choose a real number A satisfying 0 < A < — and let 
f{x,u,d{f)) ^ _ 2 

c = — — A. Then, from (14), (18) and (20) , we have 

V (x, t) <p{V (x, t)){ — XV (x, t) + fi{u)) 

< - ^p{V(x, t))V{x, t) + p{jfi{u))fi{u). y 0 < - ^Piy(x, t))V (x, t) + fi{u) 

= - (A + c)p{V (x, t))V (x, t) + fi{u) 

For this purpose, consider the following two cases. < — (AtA(x,t) + cq;(||x||)) + fi{u). 
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Thus the proof is completed. 


□ 

Remark 2.1 Lemma 2.1 and its proof can be viewed as 
an extension of the main result in [llj. If we let a{■) 
be any smooth positive definite function d(-) in Lemma 
2.1, then there exists a class ICoo function a(-) satisfying 
a(s) = O(s^) as s ^ O^, and q;(||x||) > d(x) for any 
X G M". Thus, if ai{-) satisfies lim^ sup < oo, then 

a{s) = 0(ai(s)) as s ^ O'*". Then we conelude that if 
ai(-) satisfies lim^ sup < oo, then, for any smooth 

positive definite function a{-), by Part (ii) of Lemma 2.1, 
we have 

Vix,t) < - (At/(x,t) + 05(11x11)) + ^(-a) 

< — {XV{x,t) + c&{x)) + P{u), y x,u. 


Since (3) is a switched control law, the closed loop system 
(23) is a switched nonlinear system. To analyze the sta¬ 
bility of system (23), we will resort to the multiple Lya¬ 
punov functions and average dwell time method [3,8]. 
From Theorem 4 of [3], if the closed loop system (23) 
admits multiple Lyapunov functions Up(xc), p G V, 
satisfying 


di(lkcll) < t/p(a:c) < 52 (||xc||), Vxc, 'ipG'P (24) 


-^fcp{xc,d{t)) <-XoUp{xc), Vxc, 'ipG'P (25) 

for some class /Coo functions 5i(-) and 52 ('); and some 
positive numbers Aq, then the origin of (23) is globally 
asymptotically stable for every a G Save[Td, No] with 


Td > and arbitrary Nq, where po = sup 

Xc^O 


a2(||a:e||) 

“i(ikcii) ■ 


Remark 2.2 From [10], a system of the form (7) that 
admits a function V satisfying the inequalities (8) and 
(9) is called exp-ISS, and the function V is called an exp- 
ISS Lyapunov function of (7). Moreover, by Proposition 
8 of [10] or Theorem 3 of [12], system (4) is ISS if and 
only if it is exp-ISS. In this paper, we further eall (7) 
strong exp-ISS if it admits a function V satisfying the 
inequalities (11) and (13), and call V a strong exp-ISS 
Lyapunov function of (7). Thus, Lemma 2.1 shows that 
exp-ISS is equivalent to strong exp-ISS. It will be seen in 
the proof of Theorem 3.1 that Lemma 2.1 plays the key 
role to eliminate the growth condition of the nonlinear 
functions Gi{Zi, Ci, d{t)). 


Assumption 1 For a given compact subset ID) C , 
the subsystem Zi = Fi{Zi,ei,d{t)) admits a exp-ISS 
Lyapunov funetion Vi{Zi) such that, for any d{t) G D, 

al^i\lZ^\\)<V,(Z,)<a2^i\lZ^\\), V Z, (26) 

MZr)<-X2V,iZ,)+7r,(e,), VZi,e, (27) 

where aii(-) anddi 2 i{-) are some class Koo functions with 
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aii(-) satisfying lim sup - ^, < oo, A 2 is some positive 

s^0+ ^ 

real number, and Tr^)-) is some smooth positive definite 
function. 


3 Main Result 

In this section, by combining Lemma 2.1 with the multi¬ 
ple Lyapunov functions and average dwell time method, 
we will design a distributed switched output feedback 
control law to solve the cooperative global robust stabi¬ 
lization problem for system (1). 

For convenience, let Z = col{Zi, • • • , Z^), e = col(ei, • • •, 
Bn), ei,c,(t) =col(e„o.(t)i, • • • , eycr{t)N)- From (2), we have 
Bv(T{t) = H(T{t)B for all / > 0. Then we can put the closed 
loop system composed of (1) and (3) into the following 
form 

Z = F{Z,e,d{t)) 
e = G,^t){Z,e,d{t)) 

where F{Z,e,d{t)) = col(Fi(Zi, ei, (/(/)),•• • ,F]\f{Z]\f, 
eN,d{t))), G„(^t){Z,e,d{t)) = col(Gi(Zi, ei, d(i)) + 

&ifci(e«CT(t)i), ■ • • ,GN{ZN,eN,d{t)) + bNkN{eya(t)N))■ 
Let fca(t){xc,d{t)) = co\{F{Z,e,dit)),Ga(t){Z,e,dit))) 
and Xc = co\{Z, e). Then we can further put system (22) 
into the following compact form 

Xc fca{t){Xc, d{ff) . (^I^) 


Remark 3.1 Assumption 1 implies that the subsystem 
Zi = Fi{Zi, Ci, d(t)) is exponentially input-to-state stable 
with a as the input and Vi{Zi) is an exp-ISS Lyapunov 
function of the subsystem Zi = Fi{Zi,ei,d{t)). 

Now we describe our main result as follows. 

Theorem 3.1 Under Assumption 1, for every a G 
Save[xd, Nq] with Td > and arbitrary Nq, there ex¬ 

ists a distributed switched output feedback control law of 
the form 

where pi{-), i = 1 , • • • ,N, are some sufficiently smooth 
positive functions, Xq and po are some positive real num¬ 
bers, that solves the cooperative global robust stabilization 
problem of system (1). 

Proof: First note that, since Gi{Zi,ei,d{t)) is smooth 
and Gi(0, 0, d{t)) = 0 for all d{t) G by Lemma 7.8 
of [ 6 ], there exist some smooth positive definite functions 
li{Zi), Xi(ei), such that, for all Zi G G R and 

d{t) G D, 

|G,(Z„e„d(/))|2 <+(^.) + X,(ei). (29) 
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Without loss of generality, assume bm < h < bM- By 
Lemma 2.1 and Remark 2.1, under Assumption 1, the 
subsystem Zi = Fi{Zi,ei_,d{t)) admits a strong exp- 
ISS Lyapunov function Vi{Zi) such that, for any d{t) G 
D, 

au{\\Z,\\)<V,iZ,)<a2^i\\Z,\\), V Z, (30) 

V,{Z,)<-{Xo^{Z,)+co-f^{Z,))+^^,{e,), V e„ (31) 

for some class /Coo functions dii(-) and a 2 i{-), some pos¬ 
itive real numbers Aq and cq, and some positive definite 
smooth function 

For each p £ V, let e^p = HpC and G{Z,e,d{t)) = 
co\{Gi{Zi, ei, d(t)), • • • , GAr(Zjv, cat, d(t))). Then, from 
(1), we have 


e„p = HpC = HpG{Z, e, d{t)) + HpBu 

= Gp{Z,eyp,d{t)) + HpBu ^ ^ 

where B = diag(6i, • • • , ^at), u = col(wi, • • • , un)-, and 
Gpi^Z , Cyp , d(t)) — col(Gpi i^Z , (^yp , d(t)), • * • , Gpj\[ (^Z , Cyp^ 

d{t))) = HpG{Z,H-^eyp,d{t)). Since Hp = [/iF] 
is a constant A4 matrix, then by Lemma 2.5.3 of 
[4], there exists a positive definite diagonal matrix 
Dp = diag(dpi,-- - ,dpN) such that DpHp + Hp Dp is 
positive definite. Thus B{DpHp +Hp Dp)B is also a pos¬ 
itive definite matrix. Define du = maxj max dyi\, 

peV 

dm = min{ min dpi} and Ai = {DpHp+ 

pGV ,N pGV 

HjDp)B)}. 


Let piie-ypi) = kujiie^pj), where fe is a positive real num¬ 
ber and > 1, i = 1,2, ••• ,A^, are some smooth 

non-decreasing functions to be determined later. Let 


the time derivative of 14p(e) is given by 


dVy, 


ep 


de 


iG{Z, e, d(t)) + Bu) 


dVe 


de. 


(G'p(^, Cyp, d{t)) + HpBu) 

vp 
N 


—2 ^ ^ dpib-iiOi{^Gypj^^Cypi ( Gpii^Z^ €yp^ 
N 

— ^ ^ h^jbjkLUj{eypj)ey 


N 


1 ^2 • T 

i—l i—l 

- HelpfBiDpHp + HlDp)Be:p 


N 


N 


< - ^(fcAi - £2dlib\i)ujj{elpi)e 

H- \\Gp{Z, eyp,d{t))\\‘^ 

^2 
N 


2 

vpi 


< - ^(fcAi - e2dliblj)uj‘^{elpi)e 
N 

■^|G,(Z.,e„d(t))| 


2 

vpi 


2=1 

\\H. 


pi 

£^2 


2=1 


N 


< - ^(fcAi - £2dMblj)uj‘f{elpi)e. 


2 

vpi 


2=1 


N 


+ LzlL ^('Yj(Zi) -I- Xi(ei)) 

£2 ^ 

2=1 


(36) 


Let V{Z) = y^i{Zi)- By inequality (31) we have 


Ve 


N 

i=l 


■f 


uji{s)ds. 


(33) 


N 


N 


V(Z) <-J2 (^oV^(Z,) + co7^(Z^))+J2^^((^^)■ (37) 


Then, it can be seen that t4p(e) is positive definite and 
radially unbounded. Thus, there exist some class /Coo 
functions /3ip(-) and /32p(') such that 

/3ip(||e||) < Kp(e) < /32p(||e||), Ve, Vp € V. (34) 

Choose two class /Coo functions /3i(-) and $ 2 (-) such that, 
for all p G T" and all e G R^, /3i(||e||) < /3ip(||e||) and 
P 2 [\\e\\) > /32p(||e||). Then 

/3i(||e||) < Rep(e) < ^ 2 (||e||), Ve, Vp G V. (35) 

Let e*yp = col(wi(e2pi)e„pi, • • • ,uJN{elpj^)eypN)- Then, 
by (28), (29), and (32), for any £2 > 0 and any d{t) G D, 


Finally, let Up{xc) = V{Z) -|- 14p(e). Clearly, there exist 
two class /Coo functions cii(-) and 52(-) such that the 
condition (24) is satisfied for all p G V. Also, according 
to (36) and (37), we have 

BU 

^fcp{Xy,d{t)) 

OXc 

W ^ IIH- ||2 

<-^AoC.(Z,)-5^(co-^^)7,(^.) (38) 

2=1 2=1 

N 

- “ ^ 2 dMbl^)uj'^{elpi)elpi + /5(e) 

2=1 


5 


with /3(e) = ■ By lemma 7.8 








of [ 6 ] again, for each p € V, there exist some smooth 
positive functions Pipie^pi) > 1, i = 1, ■ ■ ■ ,N such that 

p{e) = p{H~^evp) < Pip{<ivpi)elpi. Choose some 

smooth functions Pi{eypi) such that Pi^e^pi) > Pip{eypi) 
for all p € V. Then /3(e) < Pz{evpi)elp^ for all 

p GV. Define a function A/ : K+ —K such that Ai(s) = 

max Pi{eypi). Clearly, Ai(-) is non-decreasing and 
0 *^ I ^vpi I ^ ® 

Ai(|e„pj|) > Piicypi) > 1 for any e„pj S K. Then we 
choose the positive real number k such that 

k > {e2d\^b\^ + Ao^m^m + 1) (39) 

Ai 

and choose £2 such that cq — > 0. Then, 

dx d{t)) 

N N 

< - Ao ^ Vi{Z^) - Ao ^ dMbMi^‘i{elpi)elpi (40) 
2 = 1 2 = 1 

N 

-^{^Kfiipi) - Pi{e-vpi))elpi, 'ip&v 
i=l 

Finally, choose the smooth non-decreasing function uJi{-) 
such that 

(41) 

1 +e^ 

Then, since Ai(—|(£i) > Ai(|e„pj|) > Pi{eypi) > 1, we 
have 

^‘^{(iip^) > Wz(e?pJ > > Pi{eypi) (42) 


and thus 
r)U 

g/cp(Xc,d( 0 ) 

N N 

< - Ao ^ Vi{Zi) - Ao ^ dMbMi^i{elpi)elpi 

2=1 2 = 1 

< - AotA(Z) - AoKp(e) = -AoUp, Vp e V. 


Let 


po = sup 

Xc^O 



The proof is thus completed by 


invoking Theorem 4 of [3] as rephrased at the beginning 
of this section. □ 


Remark 3.2 From the proof of Theorem 3.1, we know 
that if we replaee the exp-ISS condition (27) in Assump¬ 
tion 1 by the strong exp-ISS condition (31), then the con¬ 
dition lim sup . ^, < 00 in Assumption 1 can be re- 
moved. 


Remark 3.3 The recent paper [If ] handled the coopera¬ 
tive output regulation problem for nonlinear systems with 
static directed graph. The Lyapunov function (33) here 
is similar to that used in Lemma 5 of [If]. 


Example 1 Consider the following controlled Lorenz 
multi-agent systems taken from [2[: 

Z 12 — 4“ I'li^i 

Z 2 i = L 2 %Z 2 i + ZuCi (44) 

C-i — L^2Z}_i Ci ZiiZ2i “t“ biUi, I — 1,2,3 

where bi = 1, Li = col{Lii, L 2 i, L^i) is a constant 
parameter vector that satisfies Lu > 0, L2i < 0 
and L^i > 0. Clearly, system (44) is in the form of 
(1) with Zi = col{zii,Z 2 i). To account for the un¬ 
certainty, for i =_ 1,2,3, let Li = Li -\- di where 
Li = col{Lii,L 2 i,L 3 i) = co/(3, —3.2,1.6) denotes the 
nominal value of Li, and di = col{dii,d 2 i,d^i) repre¬ 
sents the uncertainty of Li. Assume d = col{di,d 2 ,d^) S 
D = {d e \dji\ < 0.2,i,j = 1,2,3}. Let Hi = 


2 0-1 


10 0 


-1 1 0 

1 H2 = 

0 2-1 

which are both A4 

-10 2 


-1 -1 2 



matrices. Next, let a Lyapunov function candidate be 
Vi{Zi) = \z[i -\- \z\i -\- \z 2 i. Then it is possible to show 
thatVi{Z,) < -2.3zli-1.8zfi-2.bzli-^b.2e[ + 2Q.5ef < 
—A.6Vi{Zi) -I- 5.2ef -|- 26.5e^. Thus Assumption 1 is 
satisfied. 



0 2 4 6 8 10 



Fig. 1. State trajectories 


As a result, using the procedure introduced in this 
section, the cooperative robust stabilization problem 
for system (44) is solvable by a switched control law: 
u, = -12(e^^(()^ -h l)e„,^(t)i, i = 1,2,3, for any 
o-(t) € Save[Td, Nq] with Td > 2.72 and arbitrary Nq. 
Simulation is conducted for the following specific switch- 
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ing signal 


\ 1, if sT <t < {s + \)T 
a{t) =\ ^ 2 ( 45 ) 

1 ^ 2 , ^f{s+-)T<t<{s + \)T 

where s = 0,1,2, •••, and T = 2Td = 6 see. Other 
data for the simulation are di = [—0.2,0.1 ,—0 .2]^ for 
i = 1,2,3, Zi(0) = [2.6,-0.7]^, ^2(0) = [-2.8,0.9^, 
^ 3 ( 0 ) = [-1.4, ande(O) = [0.3,0.2,-0.1]'^. Figure 
1 shows the state trajectories of the closed-loop system. 
Due to the space limit, the details for designing the con¬ 
troller are omitted. 

4 An Application 

In the past few years, the cooperative control problems 
for nonlinear multi-agent systems have been extensively 
studied for the static network case in [1,2,7,13,14,15], 
and for the switching network case in [9,16]. Note that 
the nonlinear systems considered in [16] contain no dis¬ 
turbance and uncertainty, and the nonlinear systems 
considered in [9] need to satisfy certain growth condition. 
Also, the switching network is assumed to be undirected 
in [9]. In this section, we will apply our main result to 
the cooperative global robust output regulation problem 
for a class of nonlinear multi-agent systems with switch¬ 
ing network. The systems studied here contain both dis¬ 
turbance and parameter uncertainty, and do not need 
to satisfy the growth condition, moreover, the switching 
network is relaxed to be directed. 

Consider the nonlinear multi-agent systems in normal 
form with unity relative degree as follows, which is the 
same as those studied in [2,9]. 

ii =fi{zi,yi,v,w) 

yi=hi{w)uiZr gi{zi,yi,v,w) (46) 

Ci =yi - q{v,w), i = 1, • • ■ ,N 

where, for i = 1, • • • ,N, {zi,yi) G x K is the state, 
Ui G M is the input, G K is the error output, w G M"” 
is an uncertain parameter vector, and v{t) G M”” is an 
exogenous signal representing both reference input and 
disturbance. It is assumed that v{t) is generated by a 
linear system of the following form 

v = Sv, yo = q{v,w) (47) 

and all functions in (46) and (47) are globally de¬ 
fined, sufficiently smooth, and satisfy fi{0,0,0,w) = 0 , 
5 i( 0 ,0, 0, w) = 0 and q{0, w) = 0 for all w G R"”. 

As in [9], the plant (46) and the exosystem (47) to¬ 
gether can be viewed as a multi-agent system of (A -|- 1 ) 
agents with (47) as the leader and the N subsystems 


of (46) as N followers. Given the plant (46), the ex¬ 
osystem (47), and a switching signal a{t),_we can de¬ 
fine a time-varying digraph Gcr{t) = (^i^cr(t))! where 

V = {0,1 • • • , N} with 0 associated with the leader sys¬ 
tem and with i = 1, • • • ,N associate^! with the N fol¬ 
lowers, respectively, and C V x V for alH > 0. For 
all t > 0, each j = 0 , 1 , • • • , A, i = 1 , • • • , A, i 7 ^ j, 
(j, i) G £a{t) if and only if the control Uift) can make use 
of yi{t) — pj (t) for feedback control. Let Ga{t) = (V, fcr(t)) 
where V = {1 ■ • • , A}, £a{t) C V x V is obtained from 
fcr(i) by removing all edges between the node 0 and the 
nodes in V for all t > 0. 

Let A{t) = [dij{t)] G be the adjacency 

matrix of the digraph Gcr{t) i where da (t) = 0 and dij (t) = 

1 (j, i) G ^^(t); hj = 0; li ■ ■ • 1 Define the virtual 

regulated output as follows: 

N 

^va{t)i — ^ ^ dijiifjijji Vj')' 

3=0 

Let e^cr(i) col(e.„(j(t) 1 ; *' * 5 ^va{t)N^^ ^ col(ei, * * * , e^), 
eo = 0 and = [/iy(t)]"j=i with hii{t) = 

^f=odij{t) and hijft) = —dij{t) for i A j- It can be 
verified that e^cr(t) = H(T{t)^- Then our control law will 
be of the following form 

^2 — eycr(t)i'): Vi ~ fjiidlii ^V(7{t)i) 

where the functions ki and gi are sufficiently smooth 
vanishing at the origin. A control law of the form (49) 
is called a distributed switched output feedback control 
law since ei,o.(t)i is a switching signal and depends on 
iVi ~ Hj) if only if the node j is a neighbor of the node i. 
Then we describe our problem as follows: 

Problem 2 Given the multi-agent system (46), the ex¬ 
osystem ( 47 ), a group of digraphs Gp = (y,£p) with 

V = {0,1 • • • , A}, p = 1, • • • , no, and some compact sub¬ 
sets V C R"” and W C R"” with 0 G W and 0 G V, find 
TdyNo, and a control law of the form (49) such that, for 
any v(t) G V, w G W, and any a G Save[Td, Aq], the tra¬ 
jectory of the closed-loop system composed of (46) and 
(49) starting from any initial state Zi(0), t/i(0) and ? 7 i( 0 ) 
exists and is bounded for all t > 0, and limi_>oo e(t) = 0 . 

The above problem was studied recently in [9] where it 
was shown that this problem can be converted to the 
problem studied in Section 3 under the following as¬ 
sumptions. 

Assumption 2 The exosystem is neutrally stable, i.e., 
all the eigenvalues of S are semi-simple with zero real 
parts. 
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Assumption 3 For i = 1, • • • ,N, |&i(w)| > 0 for all 
w e R””. 

Assumption 4 There exist globally defined smooth 
functions Zi : R"” x R"’"— with Zi{0,w) = 0 such 
that 


dzi{v,w) 

dv 


Sv = fi{zi{y, w), q{v, w), v, w) 


(50) 


for all {v,w) G R”” X 


i = h--- ,A. 


, . 1 , dq(v, w) 

Assumptions Let Ui(v,w) = b- (--- bv — 

av 

gi{zi{v,w),q{v,w),v,w)), z = 1, • • • , A. Then Ui(z;,w) 
are polynomials in v with coefficients depending on w. 


It was further shown in [9] that the cooperative global 
robust stabilization problem of (53) was solvable by a 
control law of the form (54) under the following three 
assumptions: 

Assumption 6 For the given compact subset!]) C R"" x 
R"”, the subsystem Zi = fi(zi, Ci, d(t) admits a func¬ 
tion Vzi{zi) such that, for any d(f) G D, 


aiidl-Zill) < Vzi{zi) < a2i{\\zi\\), V Zi (56) 

Vzi{zi) <-MVziizi) + I3i{ei), 'i Zi,ei (57) 

for some class ICoa functions au{-) anda 2 i{-) withau{-) 
2 

satisfying lim sup —< oo, some positive real num- 

S-J-0+ “It is; 

ber Xi, and some smooth positive definite function 


Under Assumptions 2 to 5, there exist integers Si, i = 
1, • • • ,N, such that, for any Hurwitz matrices Mi G 
RSiXSi^ and any column vector Ai € with {Mi,Ni) 

controllable, the following linear dynamic compensator 

f]i = Mipi + NiUi, z = 1, • • • , A, (51) 

is an internal model of system (46) [9]. Moreover, there 
exist some functions 9i : vanishing at the 

origin, and some row vectors 4'^ G R^^®’ such that the 
following coordinate and input transformation on the 
internal model (51) and the plant (46) 

Zi =Zi - Zi{v, w), fji=pi- Oifu, w) - Nibf^Ci 
e-i=yi-q{v,w), zzi = zzi - 4'i?7i, z = l,---,A 

gives rise to the so-called augmented system of the plant 
(46) and the exosystem (47) as follows. 

Zi =fi{zi,ei,d{t)) 

fji =Mifji MiNibf^Ci - Nib~^gi{zi,ei, d{t)) (53) 

Bi =gi{zi, Ci, dlt)) bi^ipi 4'i AiCi -|- biUi 

where d{f) = {v,w), fi{zi,ei,d{t)) = fi{zi Zi,ei -\- 
q,v,w) - fi{zi,q,v,w) and gi{zi,ei,d(t)) = gi{zi -|- 
Zi,ei-G q,v,w) - gi{zi,q,v,w). Clearly, /i(0,0, d(t)) = 0 
and 5i(0, 0, d{t)) = 0 for any d{t) G R"" x R"”. 

By the internal model principle as can be found from [9] 
or [2], if a control law of the following form 

Ui ■> i 1 ; * * ' ; ^ (^^) 

solves the cooperative global robust stabilization prob¬ 
lem of (53), then the cooperative global robust output 
regulation of system (46) is solved by the following dis¬ 
tributed switched output feedback controller: 

(55) 


Assumption 7 For any p G V , Qp is undirected, and 
the node 0 can reach every other node of the digraph Qp. 

To introduce the last assumption, note that, since 
gi{zi,ei,d{t)) is smooth and gi{0,0,d{t)) = 0 for all 
d{t) G D, by Lemma 7.8 of [6], there exist some smooth 
positive definite functions Si{zi), li{ei) such that, for all 
Zi G R”, e, G R and d{t) G D, 

\g^{z^,ei,d{t))\'^ <Si{zi) -! li{ei). (58) 

Assumption 8 For some real number c > 0, c6i{zi) < 
V-.,{zi). 

Assumption 8 is called a growth condition on the 
nonlinear function gi{zi,ei,d{t)) which is quite restric¬ 
tive, and Assumption 7 requires the graph Ga-{t) to be 
undirected for all t > 0 which may also be restric¬ 
tive. By making use of Theorem 3.1 of this paper, 
it is possible to remove Assumption 8 and signifi¬ 
cantly relax 7. For this purpose, let Zi = col{zi,pi), 
Gi{Zi,ei,d{t)) = gi{zi,ei,d{t)) -G h'^iPi -|- and 

Fi{Zi, Ci, d{t)) = co\lfi{zi,ei, d{t)), MiPi-\-MiNibf'^a - 
Nib~^gi{zi,ei,d{t)). Then the system (53) can be put 
in exactly the same form as (1). By Theorem 3.1, As¬ 
sumption 7 can be relaxed to the following 

Assumption 9 For anyp G V , every node z = 1, • • ■ , A 
of the digraph Qp is reachable from the node 0. 

Remark 4.1 Under Assumption 9, by lemma 4 in [5], 
Hp is an A4 matrix for any p G V. 

We now further show that, by making use of Lemma 2.1 
and Remark 2.1, Assumption 8 can be removed. For this 
purpose, it suffices to show the following lemma. 

Lemma 4.1 Under Assumption 6, the subsystem Zi = 
Fi{Zi,ei,d{t)) admits a strong exp-ISS Lyapunov func¬ 
tion Vi{Zi) such that, for any d(t) G D, 

au(\\Z,\\)<mZ,)<a2^(\lZ,\\), V Z, (59) 


ni — F iPi 

Pi = Mipi -G NiUi, z = 1, • • • , A. 





V,{Z,)<-{XoMZ,)+co^,{Z,)) + 7t,{e,), (60) 

for some class ICoo functions aii(-) and a 2 i(-), some pos¬ 
itive real numbers Aq and cq, and some positive definite 
smooth function 'S’i(-). 

Proof: By Lemma 2.1 and Remark 2.1, under Assump¬ 
tion 6, the subsystem Zi_= fi{zi,ei,d(t)) admits a 
strong exp-ISS function 14- (z^) such that, for any d{t) £ 
D, 

Q:H(l|2i||) < < d2z(||2i||), V Zi (61) 

Vzi{zi)<-{XiVs,{zi)-\-ci5i{zi))-\-Pi{ei), 'i Zi.a (62) 

for some class /Coo functions Q:ii(') and 0:21 (•)> some pos¬ 
itive real numbers Ai and c \, and some positive definite 
smooth function /3i('). 

Let % = fjJPifii, Vi{Zi) = Vz, -b Vij,, and 7i(^i) = 
Siizi) -b 11%IP, where Pi is a symmetric positive definite 
matrix. Then, by Lemma 3.1 and Remark 3.3 of [9], 
Vi{Zi) satisfies both (59) and (60). □ 

It is noted that conditions (59) and (60) are the same 
as conditions (30) and (31)) of Theorem 3.1. Thus, com¬ 
bining Theorem 3.1, Remark 3.2, and Lemma 4.1 gives 
the following resnlt. 

Theorem 4.1 Under Assumptions 2-6, and 9, for ev¬ 
ery a G Savel'^d, -^o] CL'^d arbitrary Nq, 

the cooperative global robust output regulation problem 
of system (46) with the directed switching graph Gait) 'Is 
solvable by the distributed switched output feedback con¬ 
trol law of the form (55). 

5 Conclusion 

In this paper, we have established a specific changing 
supply pair technique to analyze the exponential input 
to state stability for nonlinear systems. Then, combin¬ 
ing this technique with multiple Lyapunov functions and 
the average dwell time method, we have solved the coop¬ 
erative global robust stabilization problem for a class of 
nonlinear multi-agent systems by a distributed switched 
output feedback control law. Finally, we have applied 
this result to the cooperative global robust output reg¬ 
ulation problem for nonlinear multi-agent systems in 
normal form with unity relative degree under directed 
switching network. 
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